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2+1 flavor simulations of QCD with improved staggered quarks 



O 
O 

(N 

-t— > 
O 

O 

m 

(N 



> 
On 
On 
(N 
O 

\o ■ 
O ; 

^ : 

Oh! 
I 

Oh! 



>< 



Urs M. Heller a 

American Physical Society, One Research Road, Ridge, NY 11961-9000, USA 
Oct. 5, 2006 

Abstract. The MILC collaboration has been performing realistic simulations of full QCD with 2+1 flavors 
of improved staggered quarks. Our simulations allow for controlled continuum and chiral extrapolations. 
I present results for the light pseudoscalar sector: masses and decay constants, quark masses and Gasser- 
Leutwyler low-energy constants. In addition I will present some results for heavy-light mesons, decay 
constants and semileptonic form factors, obtained in collaboration with the HPQCD and Fermilab lattice 
collaborations. Such calculations will help in the extraction of CKM matrix elements from experimental 
measurements. 

PACS. 12.38.Gc - 12.39. Fe - 12.15.Hh - 13.20.Fc 



1 Introduction 

QCD simulations with the effects of three light quark fla- 
vors fully included have been a long standing goal of lattice 
gauge theorists. The MILC collaboration has made signif- 
icant advances towards this goal, employing an improved 
staggered formalism, "asqtad" fermions This formal- 
ism, and the ever increasing computational resources avail- 
able today, have enabled simulations with the strange quark 
mass near its physical value, and the up and down quarks, 
taken to be of equal mass, as light as 1/10 the strange 
quark mass, corresponding to pion masses as low as 240 MeV 
To control the extrapolation to the continuum limit, sim- 
ulations have been done, and are still ongoing, at multiple 
lattice spacings, a ss 0.15, 0.12, 0.09 and 0.06 fm. These 
are referred to as coarser, coarse, fine and super-fine lat- 
tices, respectively. 

All the lattices generated by the MILC collaboration 
are made available to other researchers at the "NERSC 
Gauge Connection" . They have been widely used, for ex- 
ample by the LHPC and NPLQCD collaborations for the 
study of nucleon properties, form factors and scattering 
lengths, and by the Fermilab lattice, HPQCD and UKQCD 
collaborations for the study of heavy quarkonia and heavy- 
light mesons. The last three joined efforts with the MILC 
collaboration to validate these QCD simulations. They 
compared selected, "gold-plated" , i.e. well controlled, quan- 
tities computed on the lattice with their experimentally 
well known values. Agreement within errors of 1-3% was 
found EJ. 

Such a validation is, of course, valuable for any lattice 
simulation. It is, in particular, needed for simulations with 
staggered fermions, which rely on the so-called "fourth- 
root trick" to eliminate the extra species of fermions, called 



"tastes" , present with staggered fermions, and thus to gen- 
erate the correct number of sea quark flavors. A nice re- 
view of the validity or possible problems with this fourth- 
root trick can be found in [3]. 

While the MILC collaborations has made simulations 
with lighter up and down quark masses, m;, than have 
been reached by other groups, simulations at the phys- 
ical light quark mass are still too costly, even with the 
improved "asqtad" quarks. Therefore, extrapolations in 
the light quark mass - chiral extrapolations arc still 
needed to reach the physical value. Such extrapolations 
are well understood, theoretically, based on chiral pertur- 
bation theory, XPT. For results from lattice QCD, it is 
important to include a treatment of lattice effects into the 
chiral extrapolations. This holds, in particular, for sim- 
ulation results with staggered fermions, because of the 
taste symmetry breaking effects. For this purpose, XPT 
was adapted to staggered quarks in staggered chiral per- 
turbation theory, SXPT g). 

MILC's first set of accurate results for light pseudoscalarf 
based on two lattice spacings (coarse and fine, i.e. 0.12 and 
0.09 fm), one simulation strange quark mass m' s and sev- 
eral simulation light quark masses m[ , with the SXPT 
formalism for joint chiral and continuum extrapolations, 
was published in 2004 5 . To obtain quark masses, mass 
renormalization constants are needed. The values with 
one-loop Z-factors, computed by members of the HPQCD 
and UKQCD collaborations, appeared in 

Since 2004, our simulations expanded in several ways. 
On the coarse, a — 0.12 fm, lattice a second simulation 
strange quark mass m' s has been used, allowing interpo- 
lation to the physical strange quark mass. On the fine, 
a = 0.09 fm, lattice a simulation with a lighter light quark 
mass, m[ ~ 0.1m s , has been done. A coarser lattice ensem- 
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Fig. 1. The square of the pion mass divided by the sum of the 
valence quark masses as a function of the sum of the valence 
quark masses in units of n . 

ble, with a = 0.15 fm, has been added, increasing the lever 
arm for continuum extrapolations. And, finally, a super- 
fine lattice ensemble, with a = 0.06 fm, has been started. 
A run with m\ — 0.4m' s , corresponding to m x «430 MeV 
is half-finished, and lighter-mass simulations have begun. 



2 The light pseudoscalar sector 

The SVPT fitting is illustrated in Figs. □ and Dimen- 
sionful quantities are given in units of r± = 0.318(7) fm, a 
scale related to the heavy quark potential [Tj- The quark 
masses are renormalized (at one-loop) relative to those on 
the fine lattice. The data is accurate enough and at small 
enough quark masses so that the effects of chiral logs are 
evident. The red lines are the fit functions in "full contin- 
uum QCD" (valence and sea quark masses set equal) after 
extrapolation of parameters to the continuum limit. 

Preliminary numerical results for /„■ and fx obtained 
from our most recent data are 2 

f n = 128.6 ± 0.4 ± 3.0 MeV [129.5 ± 0.9 ± 3.5 MeV](l) 
f K = 155.3 ± 0.4 ± 3.1 MeV [156.6 ± 1.0 ± 3.6 MeV](2) 
fttlU = 1.208(2X1$ [1.210(4)(13)]. (3) 

Here the numbers on the left are the new values, and those 
on the right in square brackets are from Ref. 5 . In each 
case the first error is statistical, and the second systematic. 

The lattice QCD value of /k / fn can be combined with 
experimental data to extract the CKM matrix element 
V us 0. We obtain |V US | = 0.2223(l 2 j), with an accuracy 
comparable to, and a value compatible with, the latest 
(2006) PDG value, V us = 0.2257(21) [TD| . 

2 The numbers here are updated as of Lattice 2006 [Bl- 
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Fig. 2. The pion decay constant as a function of the sum of 
the valence quark masses in units of ri. 

The up, down and strange quark masses can be de- 
termined from the masses of the tt and K mesons using 
the SXPT fits. To relate these masses to the experimental 
ones requires some continuum input for isospin breaking 
and electromagnetic effects [TT]. Details can be found in 
Rcf. 5 j. Preliminary results from our current data set, 
evaluated at scale /x = 2GeV, are: 

mfs = 90(0) (5) (4) (0) MeV [76(0) (3) (7) (0) MeV] (4) 

mjf = 2.0(0) (1)(1)(1) MeV [1.7(0)(1)(2)(2) MeV] (5) 

mfs = 4.6(0) (2) (2) (1) MeV [3.9(0) (1) (4) (2) MeV] (6) 

m s /mi = 27.2(0) (4) (Q)(0) [27.4(1)(4)(0)(1)] (7) 

m u /m d = 0.42(0) (1)(0) (4) [0.43(0) (1)(0) (8)] . (8) 

Again, new results are on the left, and earlier ones from 
Refs. are in square brackets on the right. Errors are 
from statistics, simulation systematics, perturbation the- 
ory and electromagnetic effects. The main difference be- 
tween the new and old results comes from the use of a 
two- loop mass renormalization constant |12j . compared 
to the one-loop one used previously. A non-perturbative 
mass renormalization calculation is in progress. For fur- 
ther results, including Gasser-Leutwyler low-energy con- 
stants, see Ref. [8]- 

3 Heavy-light meson physics 

In a joint effort with the Fermilab lattice and HPQCD 
collaborations, we used our full QCD ensembles to study 
properties of heavy-light mesons. For the heavy c and b 
quarks we used clover fermions with the Fermilab inter- 
pretation |13| , while "asqtad" fermions where used for the 
light valence quarks. This allows to go much closer to 
the physical light quark mass than was achieved previ- 
ously. Use of SXPT, adopted to heavy- light mesons |14|. 
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Fig. 3. The ratio R q / S = Jd yjm d//d s yJmD B as function of 
m q /m a together with the SXPT fit. The black line and extrap- 
olated point are obtained after removing 0(a 2 ) effects from the 
fit, representing a continuum estimate. 
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Fig. 4. The form factor f { +\q 2 ) for D° -> K~l+u compared 
to recent measurements by Belle |20) . 

makes the necessary remaining chiral extrapolation well 
controlled, as illustrated in Fig. 03 

We found the lattice predictions [TKHTT)] 

f D , = 249 ± 3 ± 16 MeV , (9) 
f D = 201 ± 3 ± 17 MeV , (10) 
Id J Id = 1.21 ±.01 ±.04 , (11) 

compared to the later experimentally measured values from 
leptonic decays fT7lfT5] 

f D . = 283±17±7±14MeV , (12) 
f D+ = 222.6 ± 16.7±|;| MeV . (13) 

With the Fermilab lattice and HPQCD collaborations, 
we are also computing form factors for semileptonic D — > 
ir/K and B — > ir/D decays. The heavy-to- light decay am- 
plitudes are parametrized as 

(P\V\H) = U(q 2 )(p H +p P - AY + f (q 2 )A» , (14) 

where A^ = (m 2 H — m^q^/q 2 . We can compare our cal- 
culation of f { + K) {q 2 ) for D° -> K-l+v [Sj with the recent 
measurement by Belle |2()l . 



4 Conclusions 

The MILC collaboration's QCD simulations with 2 + 1 
flavors have been designed to have all errors, in particu- 
lar from chiral and continuum extrapolations, well con- 
trolled. For many quantities they have reached hitherto 
unprecedented accuracy, so that the lattice QCD calcula- 
tions are starting to have an impact on current experimen- 
tal physics programs. For example, using only our lattice 
calculation for fic/tn as a theoretical input, the CKM 
matrix V us can be determined with an accuracy compara- 
ble to the latest PDG result. A determination of the entire 
CKM matrix with lattice QCD results based on the MILC 
ensembles of gauge configurations as the only theoretical 
input has been given in Ref. |21| . 
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